Spectroscopy of ultracold atoms by periodic lattice modulations 
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We present a non-perturbative analysis of a new experimental technique for probing ultracold 
bosons in an optical lattice by periodic lattice depth modulations. This is done using the time- 
dependent density-matrix renormalization group method. We find that sharp energy absorption 
peaks are not unique to the Mott insulating phase at commensurate filling, but also exist for super- 
fluids at incommensurate filling. For strong interactions the peak structure provides an experimental 
measure of the interaction strength. Moreover, the peak height of the peaks at hco > 2U can be 
employed as a measure of the incommensurability of the system. 

PACS numbers: 



The experimental realization of Bose-Einstein conden- 
sation in weakly interacting ultracold atoms opened the 
way to numerous exciting experiments directly probing 
fundamental effects of quantum mechanics. More re- 
cently, the regime of strongly interacting atoms has be- 
come experimentally accessible by Feshbach resonances 
l| and the advent of ultracold atoms in optical lattices 
2, 3]. In a pioneering experiment, the quantum phase 
transition from a bosonic Mott insulator to a superfluid 
has been demonstrated More generally, ultracold 
atoms in optical lattices open up the possibility to simu- 
late and explore complex quantum many-body phenom- 
ena known from electronic solids, like high-temperature 
superconductivity, in a new context 0. One major ad- 
vantage of ultracold atom systems compared to their con- 
densed matter counterparts is the tunability of the sys- 
tem parameters. In particular, fast tunability in time 
made a whole area of new non-equilibrium phenomena 
experimentally accessible 0. 0- Q. The theoretical de- 
scription of these phenomena, especially beyond linear 
response, is still lacking. To probe quantum states in 
cold atoms, various new measurement techniques have 
been developed. Light-induced Bragg scattering yields 
the dynamical structure factor Q, while RF absorption 
provides information about single-particle excitations @ . 
Time-of-flight expansion gives access to the momentum 
distribution and correlations via detection of the average 
density and noise |j. Nevertheless, there is still a lack of 
measurement techniques compared to condensed matter 
setups. 

A qualitatively new way of probing the system which 
exploits fast tunability was introduced by Stoferlc et al. 
0. They determined the excitation spectrum of ultra- 
cold bosons in an optical lattice by measuring the heat- 
ing induced by a relatively strong periodical modulation 
of the lattice height (20% of the initial lattice height). 
Three distinct peaks at different frequencies were ob- 
served in the energy absorption. Up to now, these fea- 
tures, in particular the second and the third peak, are 
not well understood. Previous theoretical studies applied 



a linear response treatment 0, Ull H2I IT3I ] in analogy 
to condensed matter systems or considered weak inter- 
actions using the Gross-Pitaevskii equation How- 
ever, the relatively strong modulation and the presence 
of a trapping potential in the experiments, which im- 
plies inhomogeneous filling [TH Il6j , demand for alterna- 
tive methods to check the validity of previous approxi- 
mations. 



In this Letter we present what is to our knowledge the 
first simulation of the experiment by Stofcrle et al. tak- 
ing the full timc-dcpcndcncc into account for reasonable 
system sizes in the regime of strong and intermediate 
interaction. Hereby, we focus on one spatial dimension 
using the adaptive time-dependent density-matrix renor- 
malization group method (adaptive t-DMRG) [T7L H^ . 
It is a numerical method that allows for real-time evolu- 
tion of quantum many-body systems out of equilibrium 
with an explicitly time-dependent Hamiltonian. For the 
case of commensurate filling we compare for weak modu- 
lations our results to the analytical results of Iucci et al. 
[Tfij to show the reliability of our method. For large mod- 
ulations saturation effects occur and reduce the height of 
the absorption peaks. We show that at low temperature 
incommensurate regions with more than one particle per 
site are crucial to reproduce the multiple peak structure 
of the experimental absorption spectrum. The appear- 
ance of the peak structure is not a specific sign of the 
Mott insulating state as widely believed, but reflects the 
gap occurring for strong repulsion between the coupled 
energy bands (Hubbard bands). For strong interactions 
the peak positions can be used as an experimental mea- 
sure of the interaction strength. We also show that the 
height of the peak at frequencies frui > 2U can be em- 
ployed as a measure of the incommensurability of the 
system. Such a measure is highly relevant for the real- 
ization of quantum computing with ultracold atoms in 
optical lattices, since for current proposals the formation 
of wide regions with commensurate filling is crucial. 
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The Bose-Hubbard model [3 03 
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describes well ultracold bosons in a one-dimensional sys- 
tem subjected to an optical lattice. Here L is the number 
of sites in the chain, and bj are the creation and anni- 
hilation operators, and fij = bjbj is the number operator 
on site j. The parameters J and U are the hopping am- 
plitude and the onsite interaction strength. The third 
term models the chemical potential or an external po- 
tential, like the trapping potential. In a system with 
commensurate filling a quantum phase transition occurs 
at a finite ratio (U / J) c between a Mott insulating phase, 
in which the atoms are strongly localized at the lattice 
sites, and a superfluid phase with delocalized atoms. For 
incommensurate filling, however, it is energetically un- 
favorable to localize all the atoms and the system re- 
mains superfluid. The parameters of the Bose-Hubbard 
model are directly related to experimental quantities. For 
larg e lattice sizes an approximate formula is given by 
Hf J/E r = (4/0r)O4/£r) (3/4) exp {-2y/%jE r ) and 
U/E r = 4V2^(a s /\){V x Vl/E?)(W. Here E r is the re- 
coil energy, a s is the s-wave scattering length, and A is 
the wavelength of the laser of the optical lattice. V x de- 
notes the height of the optical lattice and V± the strongly 
confined transverse direction |24j . For a better compar- 
ison to the experiment [|J we use a s = 5.45nm of 87 Rb 
in the F = 2, nip = 2 hyperfine state and A = 825nm. 
A periodical modulation of the lattice height, here we 
use V x (t) = Vo [1 + SV cos(u;t)], translates into a periodic 
variation of the hopping coefficient J and the interaction 
coefficient U via J{t) = J[V x (t)] and U(t) = U[V x {t)]. 

V =15E r , 5V=0.01 , N=32, L=32 
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FIG. 1: Absorption rate vs. the modulation frequency u at 
a lattice depth Vo = 15-E r ; the corresponding values of the 
Bose-Hubbard model are U « 0.71E r and J « 0.0074K, i.e. 
U/J ~ 95. The lines are analytical results obtained by linear 
response treating the hopping term as a perturbation ITol . 
The symbols are the results obtained using the adaptive t- 
DMRG. Inset: Energy absorbed versus time, ui = Q.7\E r /h 
corresponds to a modulation frequency at resonance and u = 
0.5E r /h away of resonance. 



In Fig.^the dependence of the energy absorption rate 
on the modulation frequency is shown for a system which 
is initially in the Mott insulating phase with filling n = 1 
per lattice site [2j|. We found this rate a to be inde- 
pendent of the system size before saturation occurs (see 
inset Fig. QJ. It is determined by fitting the function 
/(<) = E(t = 0) + at + bi cos(cji + b 2 ) with the fitting 
parameters a, &i, and 62 to our results before saturation. 
The rate shows a two peak structure: a large peak at 
huj w U and an approximately 20 times smaller one at 
fku ?s U/2. The peak at hu; m U, in the following called U 
peak, corresponds to particle-hole excitations by a "single 
photon" process 0, floL HH If2l | . For a small modulation 
amplitude SV = 0.01 the results for the U peak obtained 
by linear response combined with perturbation theory in 
J/U (solid lines) agree very well with our full calculation. 
The peak at tiw ~ U/2, which does not occur in linear 
response, is due to "two-photon" processes where twice 
the energy quantum is needed to generate a particle- hole 
excitation |2l| . In principle absorption peaks at frequen- 
cies which are higher multiples of U can also arise. For 
commensurate filling, our calculations demonstrate that 
they are negligible |26| . 

Up to now we considered an initially Mott insulating 
state. However, if the filling is incommensurate, the sys- 
tem remains superfluid for all interaction strengths. In 
the limit of strong interaction, energy bands (Hubbard 
bands) occur separated by an energy gap of order U, but 
the system remains superfluid, since excitations in the 
highest occupied band are gapless. Therefore, as pointed 
out before [lOl ll2L |l4j , incommensurability changes the 
low (uj -c U/h) frequency spectrum. However, the pres- 
ence of a trap and the experimental resolution makes it 
difficult to observe this change in experiment. We show 
how incommensurability also imprints clear signatures in 
the high (cj ~ 2U/H) frequency spectrum. 

For a deep lattice, see Fig. [5] (a), in addition to the 
U peak a new resonance at hu> m 2U, in the following 
called 2U peak, arises [27||. The main contribution to it 
stems from the process in which one particle hops onto 
an already doubly occupied site. The process where two 
particle-hole excitations are created is negligible. There- 
fore the 2U peak is a clear indicator for regions of the 
optical lattice with incommensurate filling. For the ho- 
mogeneous system considered here we find that even at 
zero temperature the 2U peak for an incommensurate 
filling n « 1.2 is much higher than the U peak. Finite 
temperature could generate further defects in the occupa- 
tion and thereby additional contributions to the 2U peak 
as noted in ^l| for the commensurate system [2^. The 
peak height of the 2U peak could be taken as a measure 
of the 'degree' of incommensurability present in the sys- 
tem regardless of the origin of the incommensurability. 
The same finding holds for larger values of the hopping 
parameter. The peak structure becomes more compli- 
cated [cf. Fig. [21 (b)], but the peaks at ftw w 2U and 



3 



5 V=0.01 , N=38, L=32 



V = 15E r , N=38, L=32 



(a)V =15E 




0.5 1 



1.5 



2.5 



<M 1- 



15 



10 



(b) V =6 E, 




0.5 1 1.5 
hm/U 



2.5 



FIG. 2: Energy absorption rate vs. modulation frequency ui 
at incommensurate filling N — 38 particles on L — 32 sites 
for a lattice height of (a) Vq = 15E r and (b) Vo = 6E r . 
Compared to the system at commensurate filling additional 
peaks at huj w 2U and in (b) at tiui ~ 2.6U arise. Solid lines 
are guide to the eyes. 
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FIG. 3: Saturation effects for strong modulations. The ab- 
sorption rate for a weak modulation (a) SVi = 0.01 (circles) 
and a strong modulation (b) 5Vn = 0.2 (upward triangles) are 
compared. The rate for (a) is scaled by the ratio (6V2/5V1) 2 
to remove trivial scaling. Clear saturation effects in the peak 
height and width can be seen. Saturation of the integrated 
energy absorption in time: (c) is the integrated energy ab- 
sorbed up to time t m — 600h/E r per unit time (downward 
triangles) , i.e. [E(t m ) - E(0)]/t m for 8V2 = 0.2. The cho- 
sen value of 8V2 = 0.2 and t m = 600h/E r correspond to the 
experimental values. Due to saturation effects in time it is 
smaller than the energy absorption rate shown in (a). 



huj ~ 2.6U still only appear for incommcnsurably filled 
systems. The shift in position compared to the naively 
expected positions at integer multiples of U can be under- 
stood from the change in the energy structure. Whereas 
for strong interactions narrow energy bands with a gap of 
U exist, for intermediate interactions the bands broaden 
with increasing hopping. In particular, the ground state 
of the incommensurate system shifts down strongly with 
increasing hopping strength. As a result, the energy dif- 
ference towards excited states increases which causes in 
Fig. El (b) the shift of the center of the U peak to approx- 
imately 1.2U and of the 2U peak to approximately 2.6U. 
The peak around 2.1U stems from a splitting of the U 
band in processes in which a particle hole pair is created 
on singly occupied sites and doubly occupied sites. The 
ratio of the amplitudes of the peaks depends crucially 
on the degree of incommensurability in the system. We 
have verified these findings by exact diagonalization of 
Bose-Hubbard chains up to 7 sites. 

In order to use this 2U peak as a measure of the 
amount of incommensurability, one needs to determine it 
for the strongmodulations (~ 20%) used in the existing 
experiments |f|. The adaptive t-DMRG method, which 
has no problem dealing with the strong time dependent 
modulation, is thus ideally suited to tackle this question 
whereas it is delicate to use linear response for this pur- 
pose. To illustrate the difference we compare in Fig. 
the energy absorption rate for the case of a small mod- 
ulation SV± = 0.01 and a modulation with the strength 
8V2 = 0.2. The results for 8V\ are rescaled by a factor of 



{SV2/SV1) 2 to eliminate the trivial amplitude dependence 
expected from linear response. Although the structure of 
center of the peak agree qualitatively well even for these 
high modulation strength, saturation effects occur in par- 
ticular in the height of the peaks. Our results indicate 
that linear response overestimates the actual energy ab- 
sorption. For the parameters shown here a reduction by 
a factor of approximately 1.6 takes place. A bimodal 
structure seems to appear for strong modulations which 
causes a broadening of the peaks, here by a factor of ap- 
proximately 2. We also plotted the experimentally mea- 
sured quantity, namely the heating, i.e. the integrated 
absorbed energy. In Fig. [3] we show \E(t m ) — E(0)]/t m 
for a modulation time t m = 600h/E r chosen as in the ex- 
periment. The positions of the peaks agree well in the two 
spectra but the height deviates, for the parameter shown 
here approximately by a factor of 2. This is mainly due 
to saturation effects in time (cf. inset Fig. ^) which cause 
the actual integrated absorption per unit time to remain 
lower than the rate |2flj . 

In order to relate our results to present day experi- 
ments we have to investigate how much the trap affects 
the energy absorption spectrum. We find that the posi- 
tions of the peaks are robust against the presence of a 
trapping potential [3(| and that the amplitudes remain 
of the same order of magnitude. This findings can be 
seen in Fig. ^ in which we show an example of the in- 
tegrated absorbed energy at a fixed time t m versus the 
modulation frequency for the case of the presence of an 
incommensurate region in the center of the parabolic trap 
(see inset Fig. 0} 31] . Therefore as in the homogeneous 
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FIG. 4: The integrated energy absorbed up to time t m = 
50h/E r , i.e. [£(t m ) - £(0)] for SV 2 = 0.2 in a system with 
confining potential and N=65 (squares) and in a homogeneous 
system with N=46 and L=32 (circles). The inset shows the 
initial density distribution in the trapped system. 



case the occurence of a peak at Hlo « 2(7 signals the 
presence of an incommensurate region and our findings 
enables us to compare the positions in the homogeneous 
system (Fig. [2J to the experimental data . 

In the experiment for large values of the initial lat- 
tice height narrow peaks at fku as U and Tvjj 2Z7 were 
found. These peaks broaden and shift in energy if the 
lattice becomes more shallow. For intermediate interac- 
tion strength an additional peak at Hlj w 2.6(7 appeared. 
These findings for the positions of the peaks agree excel- 
lently with our results at zero temperature provided we 
assume the presence of an incommensuratcly filled region 
in the experimental system (cf. Fig. 0). 

To conclude we have demonstrated using the adap- 
tive t-DMRG method that the measurement procedure 
gives important information about the properties of 
the bosonic system, in particular about commcnsurabil- 
ity properties and the energy levels. It will therefore be 
of great interest to extend these calculations and mea- 
surements to cold fermions in optical lattices, which have 
only recently been realized experimentally |22J. 
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